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Ferromagnetism in the Hubbard model is investigated on sc, bcc, and fee lattices using a systematic 
inverse-degeneracy (l/A/") expansion which incorporates self-energy and vertex corrections such that 
spin-rotation symmetry and the Goldstone mode are explicitly preserved. First-order quantum 
corrections to magnon energies are evaluated for several cases, providing a comprehensive picture of 
the interplay of lattice, band dispersion, and interaction effects on the stability of the ferromagnetic 
state with respect to both long- and short-wavelength fluctuations. Our results support the belief 
that ferromagnetism is a generic feature of the Hubbard model at intermediate and strong coupling 
provided the DOS is sufficiently asymmetric and strongly peaked near band edge, as for fee lattice 
with finite t' . For short- wavelength modes, behavior of a characteristic energy scale ui* ~ Tc 
(magnon-DOS-peak energy) is in excellent agreement with the Tc vs. n behavior within DMFT, 
both with respect to the stable range of densities (0.20 < n < 0.85) as well as the optimal density 
n = 0.65. However, our finding of vanishing spin stiffness near optimal density highlights the role 
of long-wavelength fluctuations in further reducing the stable range of densities. 



I. INTRODUCTION 

Band ferromagnetism involves a characteristic com- 
petition between decrease in interaction energy and in- 
crease in kinetic (band) energy resulting from unequal 
spin densities in the ferromagnetic state. The Stoner 
criterion UN^ep) = 1, which signals the onset of local 
moments and highlights the strong-coupling character, 
represents the simplest (Hartree-Fock) expression of this 
competition, with the corresponding T^^ representing a 
moment-melting temperature. The physical Tc, on the 
other hand, represents a spin-disordering temperature, 
and is typically (even at RPA level) much lower than 
T^^ due to strong transverse spin fluctuations result- 
ing from the low-energy scale of collective excitations 
(magnons), which also significantly increase beyond the 
Stoner value l/7V(eF) the critical interaction strength Uc 
required for long-range magnetic order. The requirement 
of self-consistently incorporating, in the non-perturbative 
strong-coupling regime, spin and charge fluctuations in 
the magnon description, while maintaining the essential 
spin-rotation symmetry, summarizes, from the broken- 
symmetry side, the challenging nature of band ferromag- 
netism. 

Rigorous results for ferromagnetism in the Hubbard 
model have been obtained only in a few special cases. 
These include the Nagaoka theorem for a single hole in 
a half- filled band in the singular [/ — > cxo limit, ^ Lieb 
ferrimagnetisrftS, at any [/ > for half-filled bipartite lat- 
tices having asymmetry in the number of sites per sub- 
lattice, flat-band ferromagnetisro^i^i^ in specially dec- 
orated lattices having dispersionless spectrum and spe- 
cial band filling and its extension to nearly flat-band 
systems.7-8'9'1"-" 

Recent reports of ferromagnetism in the single- 
band Hubbard model on fee lattice with next-nearest- 
neighbour hopping t' in the intermediate coupling regime. 



and over a wide range of electronic densities ; i"'^^'"'^ us- 
ing a variety of techniques such as dynamical mean 
field theory (DMFT), fluctuation-exchange (FLEX) ap- 
proximation, and modified random phase approxima- 
tion (MRPA), have renewed interest in the classic prob- 
lem of metallic ferromagnetism, and have highlighted 
the importance of lattice structure and band disper- 
sion in the origin of ferromagnetism. There is now 
substantial evidence that non-bipartite lattices such as 
the fee lattice, having asymmetric density of states 
(DOS), with Fermi energy located in the region of large 
spectral weight near the band-edge, are favorable for 
ferromagnetismffi^>i^a6^a8>i9^^^ The asymme- 

try in DOS is further enhanced by the next-nearest- 
neighbor hopping term (f), providing optimal condi- 
tion for ferromagnetism in the weak and intermedi- 
ate coupling regimes f^^ii^ Similar enhancement of fer- 
romagnetism in the intermediate coupling regime due 
to f'-induced DOS asymmetry has been reported re- 
cently for the one-dimensional lattice^^ and for the square 
lattice^i^i^i22i2^i^i^i2iiH Bipartite lattices with sym- 
metric DOS such as sc and bcc lattices are, on the 
other hand, controversial candidates in which ferromag- 
netism exists only in the strong-coupling regime, if at 
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The recent investigations of existence of ferromag- 
netism on the fee lattice have also highlighted the subtle 
dependence on the distribution of spectral weight in the 
DOS induced by t'. Ferromagnetism was found to be 
absent for = in both DMFTiS and FLEXlS stud- 
ies, although at sufficiently large interaction strength 
{U/W > 2.5) a stable ferromagnetic state was obtained 
in the MRPA study. While a stable ferromagnetic state 
was obtained in a wide density range for t' = t/A in the 
DMFT study for U above a critical value Uc ^ W m 
the intermediate coupling regime, in the FLEX study for 
t' = t/2 in the weak coupling limit {U « W), fer- 
romagnetic fluctuations were found to be strongest at 
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n = 0.2. Also, the optimal density as well as the range 
of densities over which ferromagnetism is stabilized differ 
substantially in these investigations. 

However, there are certain limitations in the ap- 
proaches (DMFT, FLEX, MRPA) discussed above. In- 
corporating only the local (Ising) spin excitations, DMFT 
ignores long-wavelength spin fluctuations and the k- 
dependence of self energy. FLEX incorporates k, uj- 
dependence of self energy, but ignores vertex corrections 
of the same order, thereby breaking the spin-rotation 
symmetry. Both DMFT and FLEX are therefore not 
in accordance with the Mermin- Wagner theorem, which 
rules out long-range magnetic order in one and two di- 
mensions at any finite temperature. While an improved 
form of self-energy correction has been incorporated in 
MRPA, the momentum and frequency dependence of ver- 
tex corrections have been ignored. 

In the context of self-energy corrections in a band fer- 
romagnet, the importance of vertex corrections in restor- 
ing the spin-rotation symmetry and Goldstone mode has 
been long recognized at a formal level in terms of a 
Ward identity explicitly connecting vertex corrections to 
self-energy corrections)^ Although the spin stiffness was 
shown to be reduced from its RPA value, spin-wave ex- 
citations for arbitrary wave- vector were not considered, 
and also contribution of Stoner excitations were ignored, 
which can have significantly adverse effect on the stabil- 
ity of the ferromagnetic state, as shown in Sections III 
and IV. 

The objective of the present paper is to quantita- 
tively investigate correlation effects in a band ferromag- 
net within a spin-rotationally-symmetric scheme compat- 
ible with the Mermin- Wagner theorem. Such a scheme 
is provided by the inverse-degeneracy (l/A/") expan- 
sion within the generalized TV-orbital Hubbard modelf^S 
where quantum corrections beyond the RPA can be 
systematically incorporated while explicitly preserving 
spin-rotation symmetry and hence the Goldstone mode 
order-by-order. This approach has been applied ear- 
lier to examine quantum corrections in the antiferromag- 
netic state of the Hubbard model,'^^ where it was shown 
that in the strong-coupling limit the 1/J\f expansion be- 
comes equivalent to the 1/2S' expansion for the spin-^ 
quantum Heisenberg antiferromagnet. Recently, first- 
order quantum corrections in the ferromagnetic state 
have been obtained diagrammatically by including all 
0{1/M) self-energy and vertex corrections^^^ Within this 
approach, the Goldstone-mode spectrum for the ferro- 
magnet results in diverging transverse spin fluctuations 
and hence absence of long-range magnetic order at any 
finite temperature in one and two dimensions, exponen- 
tially large spin-correlation length in two dimensions, 
and the usual T'^/'^ Bloch-law behavior of magnetization 
at low temperature in three dimensions due to thermal 
excitation of low-energy, long-wavelegnth modes. This 
spin rotationally-symmetric approach is thus especially 
suitable for low-dimensional, low-temperature, and long- 
wavelength studies, and provides a complimentary ap- 



proach to DMFT. 

In this paper we examine the first-order (l/A/") quan- 
tum corrections to the collective spin-wave excitations 
(magnons) in the saturated ferromagnetic ground state 
of the t-t' Hubbard model. Our investigations for dif- 
ferent lattices (sc, bcc, fee), next-nearest-neighbor hop- 
pings (i'), interaction strengths (J7), and wave-vectors 
(q) provide a comprehensive picture of the interplay of 
lattice, band dispersion, and interaction effects on the 
stability of the ferromagnetic state with respect to both 
long- and short-wavelength fluctuation modes. The orga- 
nization of this paper is as follows. In Sec. II we briefly 
review the spin-rotationally-symmetric approach to the 
transverse spin fluctuation propagator. Section HI A pro- 
vides a review of the RPA-level (classical) contributions 
(delocalization and exchange) to the spin stiffness. The 
classical exchange contribution is shown to be relevant 
in determining the magnitude of the 0{1 /J\f) quantum 
corrections, which are derived and examined for different 
cases in section HI B, evaluated separately with and with- 
out the contribution of the Stoner excitations. Results 
for large-momentum, zone-boundary modes with energy 
near the magnon DOS peak (dominant mode), which is of 
interest as this energy scale provides a quantitative mea- 
sure of Tc within the spin-fluctuation theory, is discussed 
in Section IV. Conclusions are presented in Section V. 



II. TRANSVERSE SPIN FLUCTUATIONS 

Transverse spin fluctuations are gapless, low-energy 
excitations in the broken-symmetry state of magnetic 
systems possessing continuous spin-rotation symmetry. 
Therefore, at low temperatures they play an important 
role in diverse macroscopic properties such as existence 
of long-range magnetic order, temperature dependence of 
magnetization, transition temperature, spin correlations 
etc. 

We consider the time-ordered, transverse spin- 
fluctuation propagator in the broken-symmetry state, 
which describes both collective spin-wave and particle- 
hole Stoner excitations, and is given by 



Xl,-^{t-t')^i{^o\T[Sr{t)Sj{t'mo) 



(1) 



in terms of the fermion spin-lowering and -raising opera- 
tors = ^']^((T^/2)^'i. The spin-fluctuation propagator 
in q, uo space can be expressed as 



X 



l-C/0(q,c^) 



(2) 



in terms of the exact irreducible particle-hole propaga- 
tor 0(q, uj), which incorporates all self-energy and vertex 
corrections. The inverse-degeneracy expansior^i^^ 



/,(2) 



(3) 
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FIG. 1: DOS for bcc (upper) and fee (lower) lattiees, show- 
ing the significant enhaneement in DOS asymmetry resulting 
from redistribution of spectral weight due to the next-nearest- 
neighbor hopping t' . 



systematizes the diagrams in powers of the expansion pa- 
rameter 1/7V which, in analogy with 1/S for quantum 
spin systems, plays the role of h. As only the "clas- 
sical" term 0^°^ survives in the M ^ oo limit, where 
(^("'(q, w) = x°(qj'^) represents the bare particle- hole 
propagator, the usual RPA ladder series 'xp{ci,uo)/l — 
?7x''(q, (jj), amounts to a classical-level description of 
non-interacting spin-fluctuation modes. The first-order 
quantum corrections <j)'^^^ , involving self-energy and ver- 
tex corrections at 0{\/M) level, have been recently ob- 
tained for a saturated ferromagnet^ and physically in- 
corporate such effects as minority-spin spectral-weight 
transfer, correlation-induced exchange correction, cou- 
pling of spin and charge fluctuations, and magnon damp- 
ing. 

As collective spin-wave excitations are represented by 
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FIG. 2: The RPA (classical) spin stifi'ness D° for the sc lat- 
tice at U /W = 1.5, along with the hopping and exchange 
contributions. Inset shows DOS for the sc lattice. 



poles in (2), spin-rotation symmetry requires that (p = 
1/U for q,Lu = 0, corresponding to the Goldstone mode. 
Since the zeroth-order term cp'^^^ already yields exactly 
l/U for q,Lo = 0, the sum of the remaining terms 
must exactly vanish in order to preserve the Goldstone 
mode. For this cancellation to hold for arbitrary TV, each 
higher-order term in the expansion (3) must indi- 
vidually vanish, implying that spin-rotation symmetry 
is preserved order-by-order, as expected from the spin- 
rotationally invariant form [U /7V)Sj.Sj of the interaction 
term in the generalized Hubbard model. '^'^ This cancella- 
tion has been shown explicitly at the 0{1/Af) level,— not 
only for uj = but also for all to, indicating no spin- wave 
amplitude renormalization, as expected for the saturated 
ferromagnet in which there are no quantum corrections 
to magnetization. 



III. SPIN STIFFNESS 

The spin stiffness D — u>q/q'^ in the ferromagnetic 
state, defined in terms of the magnon energy ujq for small- 
q modes, provides a quantitative measure of the stability 
of the ferromagnetic state against long- wavelength fluctu- 
ations, with negative D signalling loss of long-range mag- 
netic order. We first review the different contributions to 
spin stiffness at the RPA (classical) level as their behav- 
ior and interplay provide insight into the magnitude of 
the first-order quantum corrections and of ferromagnetic- 
state stability, as discussed in the next subsection. 
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FIG. 3: Hopping (upper left) and exchange (upper right) con- 
tributions to the RPA spin stilTness (lower) for the bcc lattice 
at U/W — 1.5, showing enhancement in stiffness and reduc- 
tion in optimal density with increasing t' . 



FIG. 4: Hopping (upper left) and exchange (upper right) con- 
tributions to the RPA spin stiffness (lower) for the fee lattice 
at U/W = 1.0, showing strong reduction in the exchange con- 
tribution in the low-density range for t' = 0.5. 



A. 0(1) level (RPA): interplay of hopping and 
exchange contributions 

The bare antiparallel-spin particle-hole propagator in 
the Hartree-Fock ground state of the saturated ferromag- 
net is given by 



^(")(q,c.)^x°(q,^) 



1 



T — 

+ - »f? 



, (4) 



where normalization over number of k states is implicit 



Ek — cA are the ferromagnetic band en- 



ergies, 2A = mil is the exchange band splitting with 
magnetization m equal to the electron density n, and the 
superscript +(— ) refer to particle (hole) states above (be- 
low) the Fermi energy ep- Expanding x°(Qi'^) for small 
q, uj, with 



S = -(ck-q - Ck) = q.Veic - l-{q.V)^ek 



(5) 



one obtains 
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(6) 



which yields the RPA (classical) spin stiffness (S(")/^) 



2^^ "^^"^A~~ 



(7) 



in d dimensions. Here the angular bracket ( ) represents 
momentum summation normalized over the number of 
occupied states (^ Sk)- terms in (7) of or- 

der t and t'^/U represent hopping and exchange contri- 
butions to the spin stiffness, respectively, corresponding 
to delocalization-cnergy loss and exchange-energy gain 
upon spin twisting. As shown in section IIIB, the dom- 
inant quantum corrections involve correlation-induced 
corrections to the exchange contribution. The stability 
of the ferromagnetic state therefore involves a competi- 
tion between the hopping and exchange contributions, 
of which the latter exhibits a subtle dependence on the 
band dispersion, distribution of spectral weight in DOS, 
and the underlying lattice structure. 

We have carried out investigations for sc, bcc, and fee 
lattices, with energy dispersions 
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e^'^'^ ~ 8t cos kx cos ky cos k^ + 2t' cos 2fc^ 

e^'='= = 4t ^ cos fc^ cos k^ + 2<' ^ cos 2k^ (8) 

p<i' /J 

and electronic bandwidth W — 12t, 16t, and 16t + At' , 
respectively. We choose energy scale such that t = \. 
The numerical investigations have been carried out for 
three values of t' (0, 0.25, and 0.5) and three values of 
interaction strength {U /W— 0.5, 1.0, and 1.5) for each t' . 
The motivation behind choice of these parameters is to 
explore the subtle dependence of quantum corrections on 
lattice structure, band dispersion, and the distribution of 
DOS spectral weight in the crucial intermediate-coupling 
regime. For the fee lattice, the saturated-ferromagnet 
condition (exchange band splitting nU is greater than 
Fermi energy measured from bottom of the majority-spin 
band) is met in a wide range of electronic densities for all 
U values considered. However, for the bipartite lattices 
(sc and bcc with t' = 0), with fc^-like dispersion at band- 
edge, the Fermi energy (~ n^/^) increases more rapidly 
than n for low densities, and hence comparatively larger 
U is required to ensure the saturated-ferromagnet con- 
dition; therefore we have considered only U/W = 1.5 in 
order to cover the maximum possible range of densities. 

The limit n ^ 1 of filled majority-spin band is physi- 
cally simple and essentially lattice independent, as seen 
from Figs. 2, 3, and 4. While the hopping contribu- 
tion to spin stiffness vanishes, the exchange contribution 
approaches a finite value, thereby destabilizing the fer- 
romagnetic state. The low-density {n ^ 1) behavior of 
the exchange contribution is, however, highly sensitive to 



the distribution of spectral weight near the band bottom. 
While the hopping term slowly increases to a finite value 
with decreasing n, the exchange contribution either in- 
creases rapidly, as for the sc (Fig. 2) and bcc (Fig. 3) 
lattices, and also for the fee lattice for t' — 0.25 (Fig. 4), 
or tends to vanish, as for the fee lattice for t' = Q and 0.5 
(Fig. 4). This subtle behavior of the exchange contribu- 
tion is due to the competition between the n-dependence 
of the numerator ((Vek)^), which is highly sensitive to 
the distribution of spectral weight, and the denominator 
(2A = nil) which is independent of any lattice details. 
As seen from Figs. 3 and 4, the enhancement of spin 
stiffness with increasing t' occurs mainly due to the sig- 
nificant decrease in the exchange contribution. 

As shown below, the dominant first-order quantum 
corrections correspond to correlation-induced exchange 
contributions, and the behavior of the RPA-level ex- 
change contribution discussed above actually provides a 
qualitative idea of the magnitude of quantum corrections 
and therefore of the ferromagnetic-state stability. 

B. 0{1/J\f) quantum corrections 

Expanding the net first-order quantum corrections^ 
t/)^^-* for small q, cj in order to obtain quantum correc- 
tions to spin stiffness, we find that there is a cancellation 
not only of terms of order 1/U and terms linear in uj, 
but also of terms linear in 6 = — (ck-q — Ck), implying 
(from (5)) no quantum corrections to the delocahzation 
contribution, as expected for a ferromagnetic state. The 
surviving contributions of order q'^ can be expressed as 
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Q — irj 



(9) 



where we have set q, w = in the energy denominators as 
all terms are already explicitly second order in q. All four 
terms in (9) represent correlation-induced exchange pro- 
cesses involving minority-spin intermediate states which 
are transferred down in energy. The spectral-weight 
transfer is a correlation effect corresponding to the possi- 
bility of a site being unoccupied by a majority-spin elec- 
tron. 



Now, since 0'^^^(q) is of the form B^^^q"^, the total 
particle-hole propagator (f) = [jj — Auj — B^^\^) + B^^'' q^ , 

yielding the corrected spin stiffness as D — (i3*^°^ — 

S«)/A 

In evaluating the fJ-integral in Eq. (9), we have in- 
cluded the contribution of both the low-energy spin- 
wave excitations as well as the high-energy Stoner ex- 



6 




FIG. 5: Behavior of the corrected spin stiffness for the fee 
lattice at U/W = 1.0, evaluated separately with and without 
the contribution of Stoner excitations, shown for t' — 0.25 
(upper) and 0.5 (lower), along with the corresponding RPA 
results. 



citations. From the spectral representation of dispersion 
relation^ the il-integral is equivalent to ], 
and the imaginary part S[ ] yields finite contributions 
in the negative Q range from both collective spin-wave 
and Stoner excitations in the magnon propagator terms 
X°(Q, — Ux'^iQ, ^) etc., and also from Stoner exci- 
tations in the three antiparallel-spin particle- hole (X^k") 
terms in (9). If only the contribution from the low-energy 
spin-wave poles are included for simplicity, the integral 
in Eq. (9) can be evaluated analytically.— As Vck is odd 
in momentum, we find that the second and third terms 
in (9) give vanishingly small contributions due to partial 
cancellation, as considered earlier 

We find that quantum corrections are strong enough 
to destabilize the ferromagnetic state for the sc lattice. 



FIG. 6: The corrected spin stiffness for the fee lattice shows 
strong sensitivity to t' , with stability restricted to an inter- 
mediate range of densities for t' = 0.25 (upper) and to the 
low-density regime for t' = 0.5 (lower). The stiffness as well 
as the stable range of densities are enhanced with increas- 
ing U. The ferromagnetic state remains unstable (negative 
stiffness) at all considered values of U for t' = and also at 
U/W = 0.5 for t' = 0.25. 

bcc lattice for t' — 0, 0.25, and also for the fee lattice 
for t' — for all considered interaction strengths and 
in the whole range of electronic density. Although the 
quantum correction decreases with t' for the bcc lattice, 
as seen from the decreasing RPA exchange contribution 
(Fig. 3), it still remains larger than the 0(1) contribution 
even for t' = 0.5 except in an extremely narrow density 
range around n = 0.5, where spin stiffness was found to 
be vanishingly small. However, the ferromagnetic state 
may be stabilized at larger U as actually found for the fee 
lattice with t' — 0.25. The absence of ferromagnetic state 
on bipartite lattices in the considered range of intermedi- 
ate coupling strength is in agreement with the existing re- 
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FIG. 7: The dominant-mode energy, corresponding to modes 
with large momenta near zone boundary, suffers a strong sup- 
pression at low densities due to quantum corrections, restrict- 
ing stability to an intermediate range of densities, shown for 
the fee lattice at V jW = 1.0 for t' = 0.25 (upper) and 0.5 
(lower). 



suits obtained using variety of techniques J^ii^iii^i^ Also, 
the absence of ferromagnetism on fee lattice for = is 
in agreement with recent DMFT and FLEX studies li^ii^ 
Hereafter, the discussion refers to the fee lattice specifi- 
cally. Fig. 5 shows the substantial reduction in spin stiff- 
ness when quantum corrections are included. Results ob- 
tained by evaluating quantum corrections separately with 
and without the contribution of Stoner excitations high- 
light the effect of Stoner excitations on the stability of 
the ferromagnetic state. The dominant effect of Stoner 
excitations arises from the fourth term in (9). While 
the collective spin-wave excitations significantly reduce 
the spin stiffness from its RPA value, the contribution 
of Stoner excitations, which is particularly important in 
the low-density regime, considerably narrows the stable 



range of densities for t' = 0.25. Indeed, for t' = 0, it is the 
Stoner contribution which destabilizes the ferromagnetic 
state in the entire density range. 

For t' — 0.25, the ferromagnetic state is unstable for 
weak coupling (JJ jW — 0.5), but becomes stable in a lim- 
ited density range for intermediate couphng ([7/14^=1.0 
and 1.5), and shows (Fig. 6) an optimization near n — 
0.5. These results are qualitatively similar to the DMFT 
results,— where the ferromagnetic state was found to be 
stable for J7 > in a relatively wider density range with 
an optimal density n — 0.65. However, spatial fluctua- 
tions have been speculated to further reduce the stable 
density rangei^ For t' = 0.5, the role of both spin-wave 
and Stoner excitations becomes less significant (Fig. 5), 
which follows from the behavior of the classical exchange 
contribution (Fig. 4), resulting in stabilization of the 
ferromagnetic state even for U jW = 0.5, and significant 
enhancement in spin stiffness (Fig. 6) with decreasing n. 

For both t' = 0.25 and 0.5, increasing U results in 
enhancement of spin stiffness as well as of the stable range 
of densities (Fig. 6). This is due to both enhanced RPA 
stiffness due to reduced exchange contribution and also 
reduced Stoner contribution. 



IV. 0{l/M) QUANTUM CORRECTIONS FOR 
LARGE-MOMENTUM MODES 



We now turn to large-q (short-wavelength) modes and 
examine quantum corrections to magnon energies for q 
near the zone boundary. As the magnon DOS is strongly 
peaked near the upper end of the spectrum, and therefore 
magnon-mode summations will be dominated, especially 
in high dimensions, by contribution of modes with energy 
near the peak energy to* , we specifically consider "dom- 
inant modes" with RPA- level energy Wq k, lo* . We have 
checked for several such q and find that the quantum 
correction is quantitatively very similar. The behavior of 
the corrected magnon energy Wq with electron filling n 
is of interest as this energy scale provides a quantitative 
measure of Tc within the spin-fluctuation theory, thus al- 
lowing for a quantitative comparison with the Tc vs. n 
behavior within the DMFT which incorporates only local 
fluctuations. 

For flnite q, the magnon energy uj^ is obtained from the 
pole condition 1 — U i?e0(q, — t^^q) 0, where 0(q, w) = 
(q, w) -h 0(1) (q, uj) . The four diagrammatic contri- 
butions to the flrst-order quantum correction (j)^^^ (Qi 
for arbitrary q, lu have been given earlier. While the 
bare particle-hole propagator (f)^^\c{,Lo) remains real in 
the relevant lo range, the quantum correction ^^^^ (q, uj) 
is complex for any flnite w < due to the coupling with 
charge fluctuations. Evaluation of the U, integral is simi- 
lar as described below (9), and both collective and Stoner 
excitations are included. 

Fig. 7 shows the behavior of the dominant-mode en- 
ergy Wq with electron flUing for t' = 0.25 and 0.5. The 
nearly similar behavior shows the relatively much weaker 
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sensitivity of large-momentum modes on details of band 
dispersion as compared to small-g modes. The strong 
suppresion at low densities, a trend already present at 
the RPA level, results in an intermediate-density range 
of stability. This behavior of the dominant-mode energy 
is in excellent quantitative agreement with that of Tc ob- 
tained in the DMFT studyj^-- both with respect to the 
optimal density {n = 0.65) as well as the stable range of 
densities (0.20 < n < 0.85). Furthermore, for a realistic 
bandwidth W ^ b eV for 3d transition metals, e.g. Ni, 
the optimal magnon energy translates to approximately 
0.15 eV, which will be further reduced due to thermal 
magnon excitation at finite temperature, bringing the 
energy scale also in good agreement with the optimal 
Tc value of ~ 0.05 within DMFT. 

The substantial decrease in magnon energy at large 
momenta (Fig. 7) has been referred to as 'anoma- 
lous softening', "^^ and has been observed experimen- 
tally in ferromagnetic manganites. Anomalous soften- 
ing has also been investigated theoretically within an 
orbitally-degenerate Hubbard model- - and the ferromag- 
netic Kondo lattice model. 

Comparison of Figs. 5 and 7 shows a marked con- 
trast in the behavior of magnon energy for long- and 
short- wavelength modes for t' = 0.5, whereas the behav- 
ior is qualitatively similar for t' = 0.25. For t' = 0.5, 
while positive spin stiffness in the low-density range in- 
dicates stability with respect to long-wavelength fluctu- 
ations, the vanishing of dominant-mode energy lo* for 
short-wavelength modes signals the spontaneous onset of 
strong local fluctuations. On the other hand, while oj* 
shows a peak near n « 0.7, indicating stability with re- 
spect to local fluctuations, the vanishing spin stiffness 
near n « 0.7 indicates loss of long-range magnetic order. 
That the stability of the ferromagnetic state with respect 
to long- and short-wavelength fluctuations can differ so 
substantially is a significant feature of our approach. The 
above results show that the loss of magnetic order and 
Tc are determined by different mechanisms in the two 
density regimes — onset of strong local fluctuations in 
the low-density range and long-range fluctuations in the 
high-density range. 



V. CONCLUSIONS 

Lattice dependence of correlation effects in a band fer- 
romagnet was investigated within a systematic inverse- 
degeneracy (l/Af) expansion which allows for self-energy 
and vertex corrections in the magnon propagator to be 
incorporated in a spin-rotationally-symmetric scheme, so 
that the Goldstone mode is explicitly preserved. The 
resulting theory is therefore in accordance with Mermin- 
Wagner theorem and hence allows for long- wavelength, 
low-temperature, and low-dimensional studies, thus pro- 
viding a complementary approach to the dynamical mean 
field theory. 

The first-order (©(l/A/")) quantum corrections to 



magnon energies, evaluated separately with and without 
the contribution of Stoner excitations, were investigated 
for different lattices (sc, bcc, fee), band dispersions {f — 
0,0.25,0.5), interaction strengths {U/W = 0.5,1.0,1.5), 
and wave- vectors (g). This provided a comprehensive 
picture of the interplay of lattice, dispersion, and inter- 
action effects on the stability of the ferromagnetic state 
with respect to both long- and short-wavelength fluctua- 
tion modes. 

The present investigation supports the idea that ferro- 
magnetism is a generic feature of the Hubbard model in 
the intermediate and strong-coupling limit provided the 
DOS is sufficiently asymmetric and strongly peaked near 
a band edge (fee lattice with t'). Quantum corrections are 
strong enough to destabilize the ferromagnetic state for 
symmetric-DOS systems (bipartite lattices such as sc and 
bcc with t' = Q), as well as for weakly-asymmetric-DOS 
systems (fee lattice with t' — 0) in the whole range of 
densities and for all considered U values, in broad agree- 
ment with DMFT results. For the bcc lattice, we find 
instability of the ferromagnetic state even for finite t' (in 
the range t' < 0.5). 

Due to the similarity between the correlation-induced 
(quantum) and the RPA (classical) exchange contribu- 
tions to spin stiffness, a quantitative understanding of 
the quantum corrections, and thus of the ferromagnetic- 
state stability with respect to long-wavelength fluctua- 
tions could actually be obtained from the RPA exchange 
contribution and its competition with the delocalization 
(hopping) contribution. 

While the hopping contribution (for which the quan- 
tum correction identically vanished due to an exact can- 
cellation) was found to vanish as n — > 1 and increase 
monotonically as ti 0, both essentially lattice indepen- 
dent, the exchange contribution was found to be highly 
sensitive to band dispersion and distribution of spectral 
weight in DOS, and therefore strongly lattice dependent. 

For the fee lattice with t' = 0.25, the ferromag- 
netic state was found to be stable in an intermediate 
range of densities near n — 0.5, with both spin stiff- 
ness and stable range of densities increasing with in- 
teraction strength. Also, the density range for stabil- 
ity against long-wavelength fiuctuations was found to be 
significantly reduced as compared to that against short- 
wavelength fluctuations, in agreement with the general 
expectation when spatial fluctuations are included.-SS 

For t' = 0.5, while stiffness and stable range of den- 
sities similarly increase with U, the enhancement in 
spin stiffness with decreasing n is in sharp contrast to 
the low-density behavior of the t' = 0.25 case, which 
demonstrated the subtle dependence of quantum cor- 
rections on details of the band dispersion. For short- 
wavelength modes near zone boundary, however, the be- 
havior of the dominant-mode energy (corresponding to 
magnon-DOS peak), was found to be quite similar for 
the two cases t' = 0.25 and 0.5, with an optimization 
near n = 0.65. The behavior of the dominant-mode en- 
ergy, which provides a quantitative measure of Tc within 
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the spin-fluctuation theory for high-dimensional systems, 
was also found to be qualitatively very similar to the 
vs. n behavior within DMFT which incorporates only 
local fluctuations, suggesting that the physically impor- 
tant quantum corrections are already incorporated at the 
first-order level. 

The case of fee lattice with t' — 0.5 was found 
to exhibit contrasting behavior with respect to long- 
and short-wavelength fluctuations. The low-density 
regime is characterized by large stiffness but vanishing 
dominant-mode energy uj* ^ indicating stability against 
long-wavelength fluctuations but loss of magnetic order 
due to strong local fluctuations. On the other hand, 
the high-density regime is characterized by positive lo* 
but vanishing stiffness, indicating locally stable ferro- 
magnetic order but spontaneous onset of long- wavelength 
fluctuations. The instability of the coUinear ferromag- 



netic state in the low- and high-density regimes therefore 
involves quite different mechanisms. 

Thermal excitation of magnons and Stoner excitations 
at flnite temperature result in magnetization correction 
due to spectral-weight transfer across the Fermi energy. 
Finite-temperature magnon renormalization due to cor- 
responding additional processes in the particle-hole prop- 
agator 0(q, bj), which will yield the thermal decay of spin 
stiffness, magnon energy, and magnetization, and thereby 
determine Tc, is currently under investigation. 
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